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. Abstract 

We consider contact processes on the hierarchical group, where sites infect other sites at 
\ a rate depending on their hierarchical distance, and sites become healthy with a constant 

recovery rate. If the infection rates decay too fast as a function of the hierarchical distance, 
then we show that the critical recovery rate is zero. On the other hand, we derive sufficient 
conditions on the speed of decay of the infection rates for the process to exhibit a non- 
trivial phase transition between extinction and survival. For our sufficient conditions, we 

■ use a coupling argument that compares contact processes on the hierarchical group with 
\ freedom two with contact processes on a renormalized lattice. An interesting novelty in 

this rcnormalization argument is the use of a result due to Rogers and Pitman on Markov 
functionals. 
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1 Introduction 



1.1 Main result 

Let A be a finite or countably infinite set, called lattice, let {aii, j))i,j£A, i^^j be nonnegative 
constants, and 5 > 0. Then the contact process on A with infection rates a{i,j) and recovery 
rate 6 is the {0, Ij'^-valued Markov process X = {Xt)t>o with the following description. If 
Xt{i) = (resp. Xt{i) = 1), then we say that the site i G A is healthy (resp. infected) at 
time t > 0. An infected site i infects a healthy site j with rate a{i,j) > 0, and infected 
sites become healthy with rate 5 > 0. It can be shown (see |Lig85 Prop. 1.3.2]) that X is 



well-defined provided the infection rates are summable, in the sense that 



sup ^ a{i,j) < oo. (1.1) 
ieA, i^i 



Usually, it is convenient to assume also that ja^| < oo, where |a'l"| is defined as in (jl.ip . but 
for the reversed infection rates a^{i,j) := a{j,i). 

A contact process may be used to model the spread of an infection in a spatially ordered 
population; see |Lig99 as a general reference. A basic feature of the contact process is that 



it exhibits a phase transition between survival and extinction. Let € A be some fixed 
site, called origin. We say that a contact process on a lattice A with given infection rates 
{a{i, i^j and recovery rate 6 > survives if there is a positive probability that the 

process started with only the origin infected never recovers completely, i.e., if 



i>S() 



[Xt / Vt > 0] > 0, (1.2) 



where 6i € {0, 1}^ is defined as Si{j) := 1 if i = j and di{j) := otherwise, and G {0, 1}^ 
denotes the configuration with only healthy sites. (In typical cases, e.g. when the infection 
rates are irreducible in an appropriate sense or if the process has some translation-invariant 
structure, this definition will not depend on the choice of the origin 0.) 
For given infection rates, we let 

5c := sup {5 > : the contact process with infection rates ,^ . 

{a{i, j))ij(zf^^ ijtj and recovery rate 6 survives} 

denote the critical recovery rate. A simple monotone coupling argument shows that X survives 
for 6 < 5c and dies out for 5 > dc. 

By comparison with a subcritical branching process, it is not hard to show that 6c < 
where \a^\ is defined below (jl.ip . For a large class of lattices, it is known that moreover 6c > 0. 
For example, this is the case for nearest-neighbor processes on infinite graphs, where a{i,j) 
equals some fixed constant A > if i and j are connected by an edge and is zero otherwise, or 
if A is a finitely generated, infinite group, and the infection rates are irreducible and invariant 
under the left action of the group jSwa07l Lemma 4.18]. On groups that are not finitely 
generated, the question whether 6c > becomes more subtle. Inspired by a question that 
came up in [SwaOS] . the main aim of the present paper is to give sufficient conditions for 
(5c > (resp. 5c = 0) when A is the hierarchical group. 

By definition, the hierarchical group with freedom N is the set 

Oat := |i = (zq, ^1, . . .) : ifc G {0, . . . , — 1}, f;^. 7^ for finitely many A:}, (1-4) 
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equipped with componentwise addition modulo N. We set 

\i\ := mi{k > : = Vm > A;}, (1.5) 

and call |^ — j| the hierarchical distance between two elements i,j € ^n- We will be interested 
in contact processes on f]jv whose infection rates a{i,j) are a function of the hierarchical 
distance between i and j only. Such infection rates may always be written as 

a{ij) := a\i_^\N-\'-^\ {i,j eflN, j) (1-6) 

where (afc)fc>i are nonnegative constants. The scaling with A^^l*^-?! in (jl.6p is chosen for 
calculational convenience. It is easy to check that in order for the infection rates a{i,j) to 
summable in the sense of (II. ip . we must assume that X^^^i ak < oo. 
Here is our main result: 

Theorem 1 ((Non-) triviality of the critical recovery rate) Let N ^ 2, let (ofc)fc>i 
be nonnegative constants such that X^fcLi ^k < oo, and let 6c be the critical death rate of the 
contact process on iljv with infection rates as in U.6\) . 

(a) Assume that 

oo 

liminf iV"''log(/3fc) = -oo, where /J^ := {k>l). (1.7) 

n=k 

Then 5c = 0. 

(b) In case N is a power of 2, assume that 

oo 

7V~'^ log(Q;fc) > — oo for some m > 0. (1-8) 

k=m 

Otherwise, assume that i f j.<gj] holds with N replaced by some real N' < N . Then 5c > 0. 

The special role played by powers of 2 in part (b) is entirely due to our methods of proof and 
has no real significance. In fact, we will carry out most of our calculations for the case N = 2 
and then generalize to the statement in part (b) by a comparison argument. Note that if the 
ak have the double exponential form 

afc = e-^' (A;>1), (1.9) 

then our results show that 6c > for 1 < ^ < and 5c = for 9 > N. There is a gap 
between the conditions ()1.7p and (jl.Sp . We guess that (jl.Sp is not necessary for 6c > 0, since 
this condition is violated when infinitely many of the a^'s are zero, while it seems unlikely 
that the latter should imply 5c = 0. We do not know if condition p.7p is sharp. 

1.2 Discussion and outline 
Motivation 

Population dynamical models (but not contact processes) on the hierarchical group have been 
studied before in e.g. [SF831 IDG93[ IDawOO] . The contact process on the hierarchical group 
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which is the subject of the present paper may be used to model the spread of an infection in a 
spatiahy clustered population. Taking humans as an example, we may think of a site (zq, «i , . . .) 
as an address, where io is the house number, ii the street, 12 the town, ^3 the state and so on. 
In this example, sites at hierarchical distance less or equal than 1, 2, or 3 from a given site are 
addresses that are in the same street, town, or state, respectively. In case the ak are rapidly 
decaying, our model describes a situation where infections between large 'blocks' of sites, such 
as towns or states, are rare, hence the infection has to overcome certain 'bottlenecks' in order 
to spread and survive in the long run. In this context, we note that another model exhibiting 
such bottlenecks is the one-dimensional contact process in a random environment with fixed, 
i.i.d. infection rates between neighboring sites, see e.g. |Lig92| . Another motivation to study 
contact processes on the hierarchical group is that they may potentially be used to estimate 
contact processes on other lattices from below, including long-range processes on Z (compare 

Apart from possible applications in population biology, we believe our results are interest- 
ing from a more theoretical point of view because of the way we prove Theorem [1] (b). Finding 
upper bounds on the critical recovery rate of a contact process is generally easier than finding 
lower bounds. In line with this, the proof of Theorem [1] (a) is rather simple, but part (b) 
is much more involved. In fact, as we explain below, there are only a few known techniques 
for proving survival of contact processes in 'low' dimensions, and none seems to work well in 
our setting. The technique we finally invented is in its essence a renormalization argument. 
As such, it is interesting in the more general program of finding rigorous renormalization 
techniques for interacting particle systems. 



Renormalization 

It has been recognized long ago that the hierarchical group is especially suitable for renor- 
malization arguments. There exists an extensive literature on the Ising model on hierarchical 
lattices (see, e.g., |Dys69 IBM871 IHHWOl] ). Moreover, the self-avoiding random walk on a 



hierarchical group with 'effective' dimension four is treated in [BEI92j . while linearly inter- 
acting diffusions and near-critical percolation on the hierarchical group have been considered 
in |DG93j and [DGOG], respectively. In these last two papers, in order to get rigorous results, 
the authors take a 'local mean field limit', meaning that they send the freedom of the 
hierarchical group to infinity and rescale to get nontrivial limits. 

The intuitive idea behind our proof of Theorem [1] (b) is easily explained. For given i = 
{io,k, . . .) G ^N, set 

Bi := {(j, io,iu...)€nN:je{0,...,N- 1}} (1.10) 

Then (i?j)jgnjv is a collection of blocks Bi C JIat, each Bi containing N sites at distance 1 from 
each other. We would like to consider Bi as a single site in a 'renormalized' lattice, such that 
Bi can be either infected or healthy. Indeed, if N is large and ai > 5, then it can be shown 
that there exists a 'metastable' state on Bi in which roughly a (1 — 5/ai)-fraction of the sites 
is infected, and that transitions from this metastable state to the all- healthy state are fast and 
happen rarely. Thus, as long as 5/ai is sufficiently small, we expect our 'renormalized' blocks 
Bi to behave effectively as a single site, with an effective 'renormalized' recovery rate 5 that is 
much smaller than the original 6. Iterating this procedure, we expect the system to be more 
and more stable as we move up the spatial scale, until, in the limit, we never get extinct. 
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It may well be that this intuition can be made rigorous in a precise way for a suitably 
chosen model in the local mean field limit — > cxd, in the spirit of [DG93t IDG06] . Our 
motivation, however, was to prove results for fixed N. The problem with renormalization- 
style arguments for fixed is that in this case, one is forced to give exact bounds on how 
stable the 'metastable' state on Bi is, and how fast transitions between this state and the 
all-healthy state are. Moreover, these bounds must be translated into similar bounds on a 
renormalized lattice, in a way that can be iterated. If one tries to do this in a straightforward 
manner this soon becomes very messy and technical. 

The solution we found for this problem is a technique the second author learned about from 
a talk by Tom Kurtz on the look-down construction for Fleming- Viot processes |DK961 IDK99| 
and that originates from Rogers and Pitman |RP81j . Basically, this is a technique for adding 
structure to a Markov process X, such that if in the enriched process {X,Y), one forgets the 
added structure Y, one obtains back the original process X. An interesting feature of this 
technique is that in the enriched process {X, Y), the process X is in general not an autonomous 
Markov process, i.e., the dynamics of X depend on Y. In practice, we will set up a coupling 
between a contact process X on the hierarchical group 0,2 with freedom 2, and an 'added-on' 
process Y that lives on a renormalized lattice and that is almost a contact process itself. In 
particular, Y can be stochastically estimated from below by a contact process Y, which is 
sufficient for our purposes. For a more detailed discussion of our methods, we refer the reader 
to Sections 13.11 - 13.31 below. 



Survival 

Since the direct aim of our renormalization argument is to prove survival, we conclude this 
section with a discussion of how survival is proved for contact processes on other lattices. 
Since most of the literature deals with nearest-neighbor processes on graphs, for which there 
is just a single infection rate, it has become customary to fix the recovery rate to 1, consider 
the infection rate as a variable, and prove upper bounds on the critical infection rate. By a 
trivial rescaling of time, we may instead fix the infection rate and vary the recovery rate, hence 
any upper bound on the critical infection rate in the traditional setting can be translated into 
a lower bound on the critical recovery rate in our setting. 

If A is an infinite (connected, undirected) graph, then it is always possible to embed a 
copy of Z in A, hence the problem can be reduced to proving survival of the nearest-neighbor 
contact process on Z. (It is often possible to do better than just embedding copy of Z in A, 
see iLigSSt Thm VI.4.1].) 

For the nearest-neighbor contact process on Z, we are aware of two independent proofs that 
6c > 0. If the recovery rate 6 is sufficiently small, then it is not hard to set up a comparison 
between the contact process on Z and oriented percolation on Z x Z_|_, with a percolation 
parameter p close to one. The problem can then be reduced to showing that pc < 1 for 
oriented percolation on Z x Z^-, which is known to follow from a Peierls argument (see [Dur88l 
Chapter 5]). 

An independent aproach for proving survival of the nearest-neighbor contact process on Z, 
which gives a better bound on the critical value, is the method of Holley and Liggett [HLTSj 
(see also |Lig85 Section IV. 1]). Their basic observation is that if there exists a translation 



invariant probability law on {0, 1} such that the process X started in this initial law satisfies 

f^F[3i&As.t.Xt{i) = l]l^^ 



IPpiGAs.t. Xt(z) = l]L „>0 (1.11) 
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for all finite vl C A, then by duality P[Xo(0) = 1] gives a lower bound on the survival 
probability of the process started with a single infected site. Holley and Liggett then explicitly 
construct a renewal measure that solves (jl.lip . Their method has been refined in |Lig95| , 
leading to the best rigorous lower bound on 6c available to date. 

For lattices different from Z, there exist other, independent methods for obtaining lower 
bounds on 6c- On 7?, one may use comparison with a stochastic Ising model. On Z*^ with 
d > 3, one may use comparison with certain linear systems; this method gives the sharpest 
known bounds in high dimensions. (For both these techniques, see |Lig85 Section VI. 4].) For 



processes on trees, there is a very simple lower bound on 6c resulting from a supermartingale 
argument (see |Lig99 Thm 1.4.1]). 



In general, one can say that proving survival for contact processes is easier in higher 
dimensions. In this context, returning to the hierarchical group, we mention the following 
fact. Let ^ = {(,t)t>o be a random walk on 17 jv that jumps from a point i to j with rates a{i,j) 
as in (|1.6p . with 

afe = iV-^'(2M (A;>1), (1.12) 
where (i > is some real constant. Then it can be shown that 

P°[Ct = 0] ^ f-"^/^ (Pilogt) as t^oo, (1.13) 

where (pis a positive, periodic, continuous real function and f{t) ~ g{t) means that f{t)/g{t) — > 
1. Thus, if d is an integer, then such a random walk is similar to a usual short-range random 
walk on Z,'^. (Indeed, this is more or less to how Brydges, Evans and Imbrie construct a 
hierarchical group with 'effective' dimension four in |BEI92| . while the scaling in |DG93| is 
chosen so as to mimic the critical dimension for linear systems, which is two.) Note that in 
particular, ^ is recurrent if and only if d < 2. 

These observations are relevant when we consider comparison with linear systems as a 
method to prove survival of contact processes on the hierarchical group. Indeed, since this 
technique depends on the transience of the underlying random walk, for processes with rates 
ak as in (|1.12p . it seems this technique can only work if d > 2. Note that our Theorem [1] (b) 
shows that 6c > for any d > 0, and in fact for processes with much faster decaying rates. 

If we forget about other 'high-dimensional' techniques, this leaves us with two known 
techniques for establishing lower bounds on the critical recovery rate that might be succesful 
on the hierarchical group: comparison with oriented percolation plus a Peierls argument, or 
the method of Holley and Liggett. 

It is not hard to set up a comparison between a contact process on 12 at and some form of 
oriented percolation on Qjy x Z_|_ (with a whole set of percolation parameters pk depending 
on the hierarchical distance), but this only moves the problem to proving that the latter 
percolates if the pk are sufficiently large. Since long-range infections are essential for survival, 
it is not obvious, and seems rather difficult, to define suitable contours which could then be 
counted and estimated in a Peierls argument. 

While we did not spend much time investigating oriented percolation on x Z+, we 
did spend a considerable amount of effort trying to adapt the method of Holley and Liggett. 
As explained in [Li g95| , the renewal measure of Holley and Liggett may be interpreted as 
a certain type of Gibbs measure with the property that in (jl.lip . equality holds if A is an 
interval. The difficult part of the proof is then to show that this equality for intervals extends 
to an inequality for general subsets ^ C A. For the hierarchical group ^l^, it is not hard 
to dream up a good analogue of Liggett's Gibbs measures and to show that (jl.lip may be 
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satisfied with equality for certain special sets. (In fact, we used blocks of sites within a given 
hierarchical distance of each other.) We were not able, however, to carry out the difficult step 
in the argument, which is to extend the equality in (jl.lip for special A (the blocks) to an 
inequality for general A C ^n- It may be that this method can be carried out succesfully; our 
failure to do so is no proof that it cannot be done. 

Outline 

After proving Theorem [T] (a) in Section [21 we present and prove our coupling of contact 
processes on Q2 in Section [3] below. A more detailed discussion of our coupling can be found 
in Sections 13. 1^ 13.31 while Sections I3.4H3.8I contain proofs. The proof of Theorem [T] (b) is 
given in Section SI Appendix [X] contains a simple, but rather tedious argument needed in 
Section 13.61 

2 Extinction 

2.1 Some general notation 

Fix N > 2 and let Q = Qjsr denote the hierarchical group with freedom A^. We introduce 
contact processes whose state spaces are finite analogues of 0,. For n > 1, set 

:= {i= (zo,...,in-i) :ifc e {0,...,A-1}} (2.1) 

and 

:= {(0)}, (2.2) 

where (0) denotes the empty sequence. We equip il" with componentwise addition modulo 
A'^. For m, n > 0, we define the concatenation io j ^ ri"*"*"" of elements i G and j E 0" by 

ioj:= (io, . . .,im-iJo, ■ ■ -Jn-i)- (2.3) 

Given < m < n, by definition, the m-block in 0" with index j E f]""™- is the set 

BmU) ■.= {ioj:ie n"'} {j G < m < n). (2.4) 

We define the set of spin configurations on by 

Sn ■■= {0, 1}^" ={x = (x(i)),60" : xii) G {0, 1}} (n > 0). (2.5) 

Note that is a set containing one element and therefore Sq = {0, 1}. For < m < n, 
i € Q""™-, and x € 5"^, we define Xi € Sm by 

Xi{j) := x{j o i) (i G r?"-'", j € Jl™, X e 5„, < m < n). (2.6) 

Note that Xj describes what the spin configuration x looks like on the m-block with index i. 

For i G 0", we define \i\ as in (11. Sp with \i\ := n if in-i 7^ 0. For given 5 > and 
nonnegative constants ai, . . . ,a„, we define infection rates a{i,j) on as in (jl.6p . and we 
call the contact process with these infection rates and with recovery rate 5 the (5, qi, . . . , a^)- 
contact process. 
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2.2 Extinction 

Proof of Theorem [1] (a) For n > 0, let X*^") be the (5, ai, . . . , an)-contact process, and set 

l{n) := E"^" [inf{t > : X^"^ = 0}] (n > 0), (2.7) 

where E'^'' denotes expectation with respect to the law of the process started in 6o (compare 
()1.2p and note that now denotes the origin = (0, ... ,0) G Q"). We will estimate l{n) by 
a very crude argument. By a simple rescaling of time, we may assume that the constant \a\ 
in (11. ip satisfies \a\ = 1. By an obvious coupling, it follows that X^"'^ may be stochastically 
bounded from above by a process in Sn = {0, 1}^" where sites jump independently of 
each other from to 1 with rate 1 and from 1 to with rate 6. Obviously, the process X^^^ 
has a unique equilibrium law, which is of product form, and if X^ denotes a random variable 
distributed according to this law, then 

nX^^^ = 0] = [-^) . (2.8) 

On the other hand, since the Markov process stays on average a time (A^") ^ in the state 
every time it gets there, one has 

F[Xi^) = 0] = ^ = (2.9) 

/(n)+7V-" l + NH{n) 

where 

l{n) := E"^" [inf{t > : X^^"^ = 0}] (n > 0). (2.10) 
Solving l{n) from (j2.8p and (j2.9p and comparing with l{n), we find that 

l{n) < l{n) = A^~"((l + <5-i)^" - l) < iV-"(l + 5-^)^". (2.11) 

Now consider our original contact process on the (infinite) hierarchical group $7 n ■ We may 
stochastically estimate this process from above by a process where infections over a hierarchical 
distance > n yield infections of a new type, in such a way that infections of different types do 
not interact with each other (in particular, sites may be infected with infections of more than 
one type). Thus, in our new process, each type evolves as a (5, ai, . . . , a„)-contact process 
in some n-block, and in addition, for each k > n, each site that is infected with this type 
establishes with rate a^N^^ {N^ — N^~^) another type at a uniformly chosen site in some 
uniformly chosen n-block at hierarchical distance k. Since at any point in time there are at 
most N"^ infected sites of a given type, and each type exists for an expected time of length 
/(n), it follows that the expected number of new types created by a type during its lifetime is 
bounded from above by 

oo 

N^l{n){l-N-^) (2-12) 

k=n+l 

In view of ()2.1ip and the definition of we may estimate this quantity from above by 

(l-iV-i)(l + 5-1)^" (2.13) 

If, for some n > 1, this quantity is less than 1, then types create new types according to a 
subcritical branching process, hence a.s. at most finitely many types are created at all time. 
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hence our contact process dies out. Taking logarithms and dividing by N^, we see that for all 
6 > there exists an n > 1 such that the quantity in (j2.13p is less than one, provided that 



liminfiV-"log((l-iV~^)(l + 5~^)^"/3n+i) <0 V(5 > 0. 

n— >oo ^ ' 

This is equivalent to 

log(l + 5"^) +liminfiV-"log(/3„+i) < VJ > 0, 

ra— >oo 

which is in turn equivalent to (|1.7p . 



(2.14) 



(2.15) 



3 Coupling of contact processes 
3.1 A coupling 

Throughout this section, we fix = 2 and consider finite {6, ai, 
il" as defined in Section [2.11 We will prove the following result. 



, a„)-contact processes on 



Proposition 2 (Coupling of contact process) Let n > I, 6 > 0, and ai, . . . , a„ > 0. Let 

X = (Xt)t>o be the (5, qi, . . . , an)-contact process started in any initial law. Set 5' := 2^5 and 



a'f, := \ak+i (/c = 1, . . . , n — 1), where ^ = f{ai/5) and f denotes the function 



f{r) := 7 - \H 



with 7:=i(3+ir) (r > 0). 



(3.1) 



Then X can he coupled to a process {Y ,Y) such that {Xt^Yt)t>{) is a Markov process, (lt)t>o 
is a {6\ a'^, . . . , a'^^^^)- contact process, Iq = Yq, Yf >Yt for all t > 0, and 



¥[Yt = y I {Xs)o<s<t] = P{Xt,y) a.s. (t > 0, y G S„_i), 
where P is the probability kernel from 5„ to Sn-i defined by (recall 112. 6\) ) 
P{x,y) ■■= p{xi,yii)) {x e Sn, y £ Sn-i), 



(3.2) 



(3.3) 



with 



\ 



/p(00,0) p(00,l 

p(01,0) p(01,l 

p(10,0) p(10,l 

V p(ii,o) p(ii,i) / 



/ 1 



\ 

1 / 



(3.4) 



The coupling in Proposition [2] achieves the intuitive aim explained in Section [1.21 namely, 
to view blocks, consisting of two sites at distance one from each other, as single sites in a 
'renormalized' lattice, which can either be infected or healthy. Indeed, (j3.2p says that the 
conditional law of Yt given Xt has the following description. First, we group the sites of Xt 
into blocks, each consisting of two sites at distance one from each other. Then, independently 
for each block, if the configuration in such a block is 00 (resp. 11), then we let the corresponding 
single site in Yt be healthy (resp. infected), while if the configuration is 01 or 10, then we let 
the corresponding site in Yt be healthy with probability ^ and infected with probability 1 — ^. 
This stochastic rule is demonstrated in Figured! The transition there has probability ^(1 — ^ 
and sites in 17^ and 17^ are depicted as leaves of a binary tree. 
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10 1 

1 1 1 1 

Figure 1: Coupling of (left) and Yt (right). The conditional probability of the transition 
depicted here is ^(1 — ^). 



It is interesting that a stochastic rule for deciding whether a block is healthy or infected 
seems to work better than a deterministic rule. We will choose the function p( ■ , 1) in (j3.4p in 
such a way that this is the leading eigenfunction of a one-level (5, ai)-contact process. Thus, 
our methods combine some elements of spectral analysis with probabilistic coupling tools. 

The next lemma (which is proved in Section 13.51 below) lists some elementary properties 
of the function / defined in (jS.ip . 

Lemma 3 (The function /) The function f defined in i3.1\) is decreasing on [0, oo) and 
satisfies /(O) = | and 

f{r) =2r-^ +0{r-^) as r ^ oo. (3.5) 
3.2 Markov processes with added structure 

In Proposition [2l the coupling between the processes X and y is of a special kind. There 
exist general results that tell us how to construct processes with conditional probabilities as in 
()3.2p . such that in addition {Xt)t>o^ on its own, is a Markov process. In the present section, we 
formulate one such result, which will then be used to construct the coupling in Proposition [2l 
Let S, S' be finite sets and set S := S x S'. Let {X, Y) = {Xt,Yt)t>o be a Markov process 
with state space S and generator G. For each y ^ S' (resp. a; € S"), we define an operator 
Gy-.M.^ (resp. G'^ : M^' ^ M^') by 

Gyfix) := Gf{x,y) where J{x,y) := f{x) {x e S, y G S' , f G R^), 
G'Jiy) := G7{x,y) where J{x,y) := f{y) (x G 5, y G 5', / G M^'). 

We say that X evolves according to the generator Gy while Y = y (resp. Y evolves according 
to the generator G'^ while X = x). In particular, if Gy does not depend on y, i.e., if Gy = G 
{y G S") for some operator G : M"^ — > M"^, then we say that X is an autonomous Markov 
process with generator G. This is equivalent to the statement that for every initial law of 
the joint process {X,Y), the process X, on its own, is the Markov process with generator 
G. The next proposition (which will be proved in Section [3.4p demonstrates that even when 
X is not autonomous, it may happen that there exists an operator G such that for certain 
special initial laws of the joint process {X, Y), the process X, on its own, is the Markov process 
with generator G. It seems that Rogers and Pitman |RP81j were the first who noticed this 
phenomenon. We will prove the proposition below by elaborating on their result. 
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Proposition 4 (Markov process with added structure) Let X he a continuous-time 
Markov process with finite state space S and generator G. Let S' he a finite set, let P be 
a prohability kernel from S to S' , and let {G'^)x£S be a collection of generators of S' -valued 
Markov processes. Define an operator G : M*^ — > M'^^'^ by 

Gfix,y) := G'Jiy) (x G 5, y G 5', / e M^'), (3.7) 

and define P : M^' ^ and P : M'^^'^' hy 

Pf{x) := P{x^y)f{y) ^rid Pf{x) := P{x,y)f{x,y). (3.8) 

y&S' y&S' 

Assume that 

GPf = PGf ifeR^'). (3.9) 

Then X can be coupled to an S' -valued process Y such that {X,Y) = {Xt,Yt)t>o is a Markov 
process with state space S x S' , the process Y evolves according to the generator G'^ while 
X = X, and 

¥[Yt = y\ {Xs)o<s<t] = P{Xt, y) a.s. (t > 0, y G S'). (3.10) 

Remark 1 If X and Y are coupled as in Proposition U then it is typically not the case that 
X is an autonomous Markov process. Nevertheless, the joint Markov process (X, Y) has the 
property that if the initial law satisfies 

P[yo = y|^o] =mo,y) a.s. (t>0, yG^'), (3.11) 

then X, on its own, is the Markov process with generator G, and (j3.10p holds. 

Remark 2 If X and Y are coupled as in Proposition HJ then it may happen that Y is an 
autonomous Markov process. In this case, we will say that Y is an averaged Markov process 
associated with X. In the general case, we will say that Y is an added-on process. 



3.3 Discussion 

We mention a few open problems concerning our coupling. 

1° Can one modify Proposition [2] such that Y = Y, i.e., (in terminology invented in the 
previous section), for a given (5, ai, . . . , a„)-contact process X, can we find a {5' , a'l, . . . , a'^_i)- 
contact process Y such that Y is an averaged Markov process of X? This would probably 
involve a kernel P and constants S',a'i, . . . that are more difficult to describe and less 

explicit than the ones in Proposition [2] but would have great theoretical value, since the 
resulting map {S,ai, . . . , an) i— > {6' , a'l, . . . , a^„^) would represent a rigorous renormalization 
transformation. 

2° Is it possible to construct a similar coupling as in Proposition [2l but with Yt < Yt? This 
could potentially be used to relax condition (|1.7p . 

3° Can one use Proposition [2] to construct a probability law on {0, 1}^^ that satisfies 
condition (jl.lip of Holley and Liggett? This would not add much in the line of proving 
survival (which is already achieved) but might add to our understanding of the method of 



Holley and Liggett, which is rather poor. In particular, in Lig95 it is shown that this method 



may be used to calculate a sequence of approximations of the critical recovery rate, but beyond 
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the second member of that sequence, there is no proof that these approximations are lower 
bounds on 5c (though they are conjectured to be so). 

4° Is it possible to make the methods of the present paper work on Z instead of 1^2? At 
first sight, it seems that the hierarchical structure of is essential to Proposition [2l However, 
when we think of the latter as 'forgetting the fast modes of the spectrum', something may be 
possible. Any link between Proposition [2] and the method of Holley and Liggett might also 
provide a clue. 



3.4 Added-on processes 

Proof of Proposition [4] We adopt the convention that sums over x,x',x" always run over 
S and sums over y,y', y" always run over S". Write 

Gf{x) = Y,r{x,x')[f{x')-f{x)), 

G'j{y) = jZ''Ay.y'){f{y')-f{y)). ^^'^'^^ 

y' 

where r(x, x') (resp. r!^(?/, y')) denotes the rate at which the Markov process with generator G 
(resp. G'^) jumps from a state x to a state x' (resp. from y to y'). 

Set 5 := S" X S' . We let {X,Y) = {Xt^Yt)t>o be the Markov process in S started in an 
initial law satisfying (j3.1ip . with generator G defined by 

Gf{x, y) := ^ ty{x, x') {fix, y) - f{x, y)) 

+ E r',{y,y'){f{x,y')-f{x,y)) 

y':P{x,y')>0 

+ J2 ''-(y' y') E 'iy (^(^'' y'^ - y^) 

y':P{x,y')=0 x' 

i{x,y) e S, f e R^), where 

w /N r{x,x')P{x',y) , r{x,x')P{x',y) 

ty{x,x):= — ^ and qy{x,x):=— 3.14 

P{x,y) l^^„r{x,x")P{x" ,y) 

These formulas are not defined if P{x,y) = resp. X^^,;; r{x, x")P{x" ,y) = 0, so in the first 
case we define ty{x,x'), in some arbitrary way, while in the second case we choose for qy'{x, ■ ) 
some arbitrary probability distribution on S. In any case, it will be true that 

{^^r{x,x")P{x" ,y)^qy{x,x') = r{x, x')P{x' ,y), (3.15) 

x" 

since the right-hand side of this equation is zero if J2x" '^(^j y) = 0- 

Formula (j3.13p says that the process {X,Y) jumps from a state {x,y) to a state {x',y) 
with rate ty{x, x'). In addition, while X is in the state x, the process Y jumps from the state 
y to the state y' with rate r'^{y,y'). During such a jump, if P{x,y') = 0, then the process 
X does nothing but if P{x, y') = 0, then the process X jumps at the same time to a state 
x' chosen according to the probability kernel qy'{x,x'). In particular, these rules say that the 
process Y evolves according to the generator G'^ while X = x. 
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It is known that (j3.10p holds for the process {X, Y) started in any initial law satisfying 
(jS.lip . provided that 

GPf = PGf (/gM^). (3.16) 

The sufficiency of (I3.16P follows, for example, from |Kur98t Corollary 3.5], which is a rather 
technical statement about martingale problems. A much less technical version of this result 
can be found in |RP81j . 

We may rewrite G in the form 

Gf{x, y) := ^ tyix, x') {fix', y) - f{x, y)) 

x' 

y' 

+ Y r'^{y,y')Yqy'{x,x'){f{x',y')-f{x,y')) {x e S, y e R). 

y':P{x,y')=0 x' 

(3.17) 

We calculate, remembering the definition of ty{x,x'), and letting G{x,x') denote the matrix 
associated with the operator G, 

GPfix) = rix, x') ( Y P(^'^ y)f(^'^ y)-J2 ^(^' y)) 

x' y y 

= E E [^(^'' (/(^'' - y)) + (^(^'' - y)] 

y x' 

= E y^ E (■^(^'' ~ y^) (3.18) 



y. P{x,y)=0 x' 



and 



PGfix) = J2 Pix, y) E *y (/(^'' y) - /(^' y)) 
+ Y ^(^' y') (/(^' y') - /(^' y)) 

yy' 

= Pi^^ y) E (^(^'' y) - y)) 
+ E E ( E ^(^' y'y- (y'' y)) (/(^'' y) - /(^. y)) 

y:P(x,j/)=0 x' y' 

+ E(E^(^'y')Gx(y',y))/(x,y), 



(3.19) 



s/ y' 

where to get the second equality we have reordered our terms and relabelled indices. The first 
terms on the right-hand sides of (I3.18|] and ()3.19p are equal while the third terms agree by 



13 



()3.9p . Since by (j3.9p . for each x,y such that P{x,y) = 0, one has 

J/' 1^' x" x" 

(3.20) 

we see by ()3.15p that also the second terms on the right-hand sides of ()3.18p and (I3.19P agree, 
hence (j3.16p holds. I 

Remark For fixed y G S", set Sy := {x ^ S : P{x,y) > 0} and consider the operator Gy 
defined by (compare ([333]l - ([3Tl]l ) 

Gyfix) := ty{x,x')[f{x') - fix)) {x eSy, fe R^y). (3.21) 

Then Gy is a 'compensated /i-transform' of the operator G, with the function h{x) := P{x, y). 
Here, if G is the generator of a Markov process on S and h is a nonnegative function on S, 
then 

G^f ■.= h-^G{hf)-h-^{Gh)f (3.22) 

defines a generator of a Markov process on the space Sh '■= {x : h{x) > 0}. This sort of 
transformation has been called a compensated h-transform in |FS04j . In particular, if h is 
harmonic, i.e., Gh = 0, then G'* is the usual /i-transform of G. 

3.5 Definition of the added-on process 

In this section, we prove Proposition [2j Our proof depends on some calculations that will be 
done in the next three sections. We wish to construct an S'^^^-valued added-on process Y 
on X, such that Y can be stochastically estimated from below by a (6' , a\, . . . , a^_]^)-contact 
process. We introduce the notation 

x(z,j) := (x(i),x(i)) (x G Sn, i,j G n^)- (3.23) 

With this notation, the generator of X can be written as 

Gf{x)=6 lM^)=l}{f{x - 5,) - fix)) 

ien" 

n 

+ 1Mm)=(o,i)}(/(^ + '^0-/(^))- ^^-^^^ 

\i~j\=k 

For any x G 5i = {0, 1}^, we write 

'00 if X = (0,0), 
x :=< 01 if X = (0,1) or (1,0), (3.25) 
. 11 ifx = (l,l). 

For each x £ Sn, we define a generator G'^ of an S'„_i-valued Markov process by (recall (|2.6p l 

G'xf{y) = s' Y hy{^)=l}{f(y-^^)-f(y)) 

n-l 

+ ^afc+i2^'= Y Hxi,Xj)l^y(^ij)=(^o,i)} (3.26) 
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where 6' is defined as in Proposition [2] and a, b are the functions 
a(00,00) a(00,01) a(00,ll) 

a(01,00) a(01,01) a(01, 11) 1=1*^ 1 | (3-27) 





a(ll,00) a(ll,01) a(ll,ll) 

and 

6(00,00) 6(00,01) 6(00,11) 

6(01,00) 6(01,01) 6(01,11) I = I i * I . (3.28) 
6(11,00) 6(11,01) 6(11,11) 

Here ^ is defined as in Proposition [2] and the symbol * indicates that the definition of a and 
6 in these points is irrelevant. Indeed, in the next three sections, we will prove the following 
fact. 

Lemma 5 (Added-on process) Let G he the generator of the (5, ai, . . . , an)-contact process 
on Sn, let P be the probability kernel from Sn to Sn-i defined in Proposition\^ and let {G'^)x£Sn 
be the generators defined in i3.26\) . where the functions a and b are defined as in { 3. 27^ - ^3. 28\) . 
Then, no matter how we define a and 6 in points indicated with the symbol *, one has 

GPf = PGf if e M'^"-!), (3.29) 

where G,P, and P are defined as in ^3. 7\ )- f37^) . 

Based on Lemma [5l we can now prove Proposition [2j 

Proof of Proposition [2] By Proposition U] and Lemma[5l we can couple X to an S'^-i-valued 
process Y such that {Xt,Yt)t>o is a Markov process, Y evolves according to the generator G'^ 
while X = X, and ()3.2p holds. 

By Lemma [3l < ^ < ^. It follows that the functions a and 6 in (j3.27p - (|3.28p satisfy 
a > ^ and 6 > 0. From this and (j3.26p . it is easy to see that {X,Y) can be coupled to a 
{6' , a[, . . . , a^_]^)-contact process Y, such that Yq = Yq and Yt >Yt for all t >0. I 

For completeness, we give here the: 



Proof of Lemma [3] Set ^(7) := 1 ~ \J l"^ ~ 2- Then it is straightforward to check that 
= \. Moreover, |^C(7) = 1 - 7(7' " \r^'^ = 1 - (1 - < on [f,oo), so 



7 h- > ,^(7) is decreasing on [|,oo). Set e := 7 ^. Then (^{e ^) = e -"^(l — ^1 — ^e^). We 
observe that 

(1 - ^1 - i.^) 1^^^ = 0, 



■(i-y^)L^o = fHi-l-^)-^/^l 



(3.30) 



de\2^ 2^ > ^4^ 2^ ' ) \e=0 2' 

hence i{e-^) = e^^{\e'^ + 0{e^)) = \e + 0{e'^), i.e., 

e(7) = i7"' + 0(7"') as 7^00. (3.31) 

To translate this to the statements in Lemma [3l it suffices to note that the function r 1— > 
7(r) := 1(2 + |r) is increasing on [0, 00), satisfies 7(0) = |, and 7(r) = \r + 0(1) as r ^ cxD.i 
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3.6 Reduction to a one- and two- level system 

In this section, we prove Lemma [5l Our proof is based on two lemmas which wiU be proved 
in the next two sections. 

Proof of Lemma [5] We start by rewriting the generator in ()3.24p as follows: 

Gf{x) =S l{x(i'oi)=i} {fix - 6eoi) - fix)) 

+ai2^^ ^ ^ l{x{i'oi)=o, x{i"oi)=i}{fix + Si'oi) - fix)) 
i^^"-' i',i"en^ 



i'-i" =1 



{x{i' oi)=0, x{j'oj)=l} {fix + 5i'oi) - fix)) (3.32) 



k=2 . i'j'eni 

\i-j\=k-i 

n-l 

Y R^f ix)+Y^k+l2-' ^^jf(' 



jgQn-l k=l 



\i~j\=k 



where 



(/GM^"-i, i/j), (3.36) 



(ii) hjPf = Phjf 



(3.33) 



Rif{x):=5 Y '^{x{i'oi)=i}ifix - 5iioi) - fix)) 
i'e{o,i} 

+ai2^^ Y^ l{x(i'oi)=o, x{i"oi)=i} {fix + Si'oi) - fix)) , 

i',j"G{0,l} 

i'yti" 

hjfix) := 2~^ Y l{x{i'oi)=o, x{j'oj)=i} {fix + 5i'oi) - fix)) . 
i',i'e{o,i} 

Likewise, we may write the operator in (I3.26P as 

n-l 

G'Jiy)= Y ^^/(y) + E«'^+i2-'^ E (3.34) 

fc=i i,jGn"-i 

\i-i\=k 

where 

i?^/(y):=<5'l{,(,)=i}(/(2/-<5,)-/(2/)), 

I%fiy) ■= [a(^j,^j)l{j/«=o, s/(j)=i} + ^(^i'%)l{s/(i)=o, y(j)=o}] (/(y + ^i) - fiv))- 
In view of (I3.32p and ()3.34p . in order to prove ()3.29p . it suffices to show that 

(i) RiPf = PE!j ] 



(3.35) 



where Uj : R-^"-! ]^SnxS„^i jg (;iefined as Iijfix,y) := I'^jfiy). Note that since R[ does not 
depend on x, there is no need to define i?,. 

The operators Ri,R^,Iij, and lij act only on certain coordinates. In view of this, our 
problem reduces to a lower-dimensional one, and (j3.36p follows from Lemmas [6] and [7] stated 
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below. It is not difficult, but notationally cumbersome, to give a formal derivation of ()3.36p 
from Lemmas [6] and [3 For completeness, we give this derivation in Appendix [Al I 

Recall that = {0, 1}^' and = {0, 1}^" = {0, 1}. Let 5 > 0, ai > 0, and let 5' be defined 
as in Proposition [2l Let P be the probability kernel from Si to 5*0 defined in (|3.3p - (j3.4p and 
let P : M*^" — > R*^! be defined as in (13. 8p . Let R be the generator of the {6, ai)-contact process 
X on 5*1 and let R' be the generator of the (5'-contact process Y on Sq. The latter is just the 
Markov process with state space {0, 1} that jumps from 1 to with rate 5' . The next lemma 
implies that Y is an averaged Markov process associated with X, i.e., X and Y can be coupled 
such that (f3T0]) holds. 

Lemma 6 (One-level system) One has 

RPf = PR'f (/gM^o). (3.37) 

Formula (13.37P implies (13.36P (i). We next formulate a lemma that implies (j3.36p (ii). 

Let 5 > 0, ai > 0, let P be the probability kernel from 5*2 to Si defined in (|3.3p - ()3.4p . and 
let a, b be the functions in (I3.27p - (j3.28p . We define a generator / of a Markov process in S2 
and generators {I'^)x<^S2 °f Markov processes in by 

^fi^) ■=h Yl ^{^M)=o, x(i,i)=i} {fix + %,o)) - fix)) , 

jjg{o,i} (3.38) 

I'xfiy) ■= [a(^o,^i)l{j;(o,i)=(o,i)} + ^(^o,^i)l{s/(o,i)=(o,o)}] {fiy + ^0) - fiy))- 

Our next lemma says that the {I'x)x<^S2 define an added-on Markov process associated with 
the process with generator /. 

Lemma 7 (Two-level system) One has 

IPf = Plf (/GM^i), (3.39) 
where P, P are defined as in 113. 8\} and If{x,y) := I'^fiv) ix G 5*2, y G 5i). 

3.7 The one-level system 

Proof of Lemma [6] We may write (|3.37p in the matrix form 

E Rix,x')P{x',y)fiy) =Y.Y. Pix,y')R'iy',y)fiy) ix e Si, f € M^«), 

x'€Si yeSo y'eSo yeSo 

(3.40) 

which is equivalent to 

E R{x, x')p{x', y)=Yl Pi^^ y')R'iy', y) ixeSi,ye So), (3.41) 

x'&Si y'eSo 



where p is the function in (j3.4p . Here 



Y Rix,x')pix',y)=Rp{-,y){x) (x G {0, 1^, y G {0, 1}). (3.42) 

x'eSi 
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Thus, (j3.4ip says that Rp{ ■ , 0) and Rp{ ■ , 1) can be written as a hnear combination of the 
functions p{- ,0) and p{ ■ , 1). It follows that T := span{p( ■ ,0),p{ - , 1)} is an invariant sub- 
space of the operator R. Since p{ ■ , 0)+p{ ■ , 1) = 1, the space J- contains the constant function 
1. We will show that J- is in fact the span of 1 and one nontrivial eigenfunction of R. 

We start by noting that by symmetry, the space H := {f ^ M'^^ : /(0, 1) = /(1,0)} 
is invariant under R. Since (0,0) is a trap of the (5, Qi)-contact process, the space Ho '■= 
{f G Ti. : /(0,0) = 0} is also invariant under R; in fact, H is the span of Hq and the trivial 
eigenfunction 1. In view of this, we look for eigenfunctions of R in TiQ. We observe that for 
/G Wo, 

RfiO, 5{0- /(0, 1)) + i«i (/(1, 1) - /(0, 1)) 

Rfihl) J V 25(/(0,l)-/(l,l)) 

(3.43) 

-(5 + iai) lai \ / /(0,1) ^ 
26 -25 j[f{l,l) 



To find the eigenvalues, we must solve 

-{6 + ^ai 
25 



det 



)-A 



-26 -X 



(3.44) 



which gives 



{6 + ^ai + X){26 + X) = 6ai 
^ + (35 + iai)A + 2(5^ 







44> 

^ X 
^ X 



A + 



35 + iai\2 /35+1q,i\2 




-26{j±^j^-^) 



(3.45) 



where 7 := |(3 + ^^) (compare (13. 1|) ). In particular, the leading eigenvalue is A = —25^ = 
—6', where ^ and 6' are defined as in Proposition [2l To find the corresponding eigenfunction, 
we need to solve 



2<5(/(0,l)-/(l,l)) =A/(1,1) 

^ 25/(0, 1) = {26 + A)/(l, 1) = 25(1 - 0/(1, 1) 

^ /(0,1) = (1-0/(1,1), 



(3.46) 



which yields the eigenfunction 



/ /(o,o) \ 




/(o,i) 


i-e 


/(i,o) 


i-e 


V /(i,i) / 


V 1 / 



(3.47) 



Our calculations so far show that := span{l,p( • , 1)} = span{p( • , 0),p( • , 1)} is an invariant 
subspace of the operator R. It follows that there exist constants {R' {y' ,y))y,yie{o,i} such that 



y'e{0,l} 



(3.48) 
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In fact 



{Rp{-,0)Rpi-,l)) = {6'pi.,l) -6'pi.,l))= (^'(•'O)^'(-'l)) r 



^ 6' -5' 

I1GIIC6 

i?'(0,0) i?'(0,l) \ _ ( ^ 



(3.49) 



(3.50) 



which we recognize as the generator of a Markov process on {0, 1} that jumps from 1 to 
with rate 6' . I 



3.8 The two-level system 

Proof of Lemma [7] We may write (|3.39p in the matrix form 

E i{x,x')P{x',y)f{y) =Y.Y. Pix,y')i'Ay',y)fiy) g ^2, / g m^^), (3.51) 

x'€S2y<^Si y'eSiyeSi 
which is equivalent to 

^ lix, x')Pix', y)= Yl ^(^' y) ixGS2, yeSi). (3.52) 

x'£S2 y'eSi 

Here 

Y x')P{x', y) = IP{-, y){x) {x e S2, y e Si). (3.53) 

Note that Sn = {0, 1}^" = {0, 1}{0'1>" has 2^" elements, so = 2^ = 4 and jSal = 2^ = 16, 
hence {IP{ - ,y){x))x<^Si, y<^S2 is a matrix with 4 • 16 = 64 entries. Luckily, using symmetry, 
we can reduce the size of our problem quite a bit. We start by calculating 

P{x,y) = Py{xo,xi) = p{xo,ymp{xi,y{l)) (3.54) 

for xo,xi,y G {(0, 0), (0, 1), (1, 1)}. For brevity, we write 00 = (0,0), 01 = (0,0), and 11 = 
(1,1). We have 

Poo(00,00) Poo(00,01) Poo(00,ll) \ /l ^ 0\ 
Poo(01,00) Poo(01,01) Poo(01,ll) = e , (3.55) 
Poo(ll,00) Poo(ll,01) Poo(ll,ll) / Vo 0/ 

Poi(00,00) Poi(00,01) Poi(00,ll) \ /O 1-e 1 
Poi(01,00) Poi(01,01) Poi(01,ll) = e(l-0 e I, (3.56) 



and 



Poi(ll,00) Poi(ll,01) Poi(ll,ll) / \0 



Pii(00,00) Pii(00,01) Pii(00,ll) \ / 
Pii(01,00) Pii(01,01) Pii(01,ll) = (1-0' 1-e |. (3.57) 
Pii(ll,00) Pii(ll,01) Pii(ll,ll) / \0 1-e 1 
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Recall the definition of x from (j3.25p . If {X{t))t>o = {XQ{t), Xi{t))t>o is a Markov process in 
S2 = Si X Si with generator /, then {XQ{t), Xi{t))t>Q is a Markov process that jumps with 
the following rates: 

(00,01) 



(00,00) 
(01,00) 
(11,00) 

From this, we see that the functions IP( 



1 

(01,01) 

1 
2 

(11,01) 



and 



/Poo(00,00) /Poo(00,01) /Poo(00,ll) 
/Poo(01,00) /Poo(01,01) /Poo(01,ll) 
/Poo(ll,00) /Poo(ll,01) /Poo(ll,ll) 

IPoi(00,00) /Poi(00,01) /Poi(00,ll) 
IPoi(01,00) /Poi(01,01) /Poi(01,ll) 
IPoi(ll,00) /Poi(ll,01) /Poi(ll,ll) 



/Pii(00,00) /Pii(00,01) /Pii(00, 11) 
/Pii(01,00) /Pii(01,01) /Pii(01,ll) 
/Pii(ll,00) /Pii(ll,01) /Pii(ll,ll) 



(00,11) 
2 

(01', 11) 
1 

(ll','ll). 
IPy{x) are given by 

-m-e) 

2? 



(3.58) 



2 


-^e(i-e) 


(i-e)^ 




-2(1-6 



(3.59) 



(3.60) 




(3.61) 



We wish to express the functions {IP{ ■ , y))y,zs^ in the functions {P{ • , y))y)zs-i ■ Unlike in the 
previous section, the span of the functions (P( • ,y))yeSi is not invariant under the operator 
/, so we cannot express the functions (/P( • ,y))y£Si ^ ^ linear combination of the functions 
(P( • ,y))y^Si- However, we can find expressions of the form (compare (|3.52p l 



IP(-,y)(x) 



5] P{x,y')lUy',y) 

y'eSi 



{x e S2, ye Si), 



(3.62) 



where the coefficients I'^^y' , y) do not depend too strongly on x. Solutions to this problem are 
not unique. The claim of Lemma [7] is that we can choose 



/ 1^(00,00) 
4(01,00) 
1^(10,00) 



( 



\ 



7^(00,01) 
4(01,01) 
4(10,01) 

4(11,01) 

-h{xQ,xi) 






4(00, 10) 
4(01,10) 
4(10,10) 
4(11,10) 





-a{xQ,xi) 





4(00,11) \ 
4(01,11 
4(10,11 
4(11,11) J 

b{xo,xi) 








a(xo,xi) 





(3.63) 



where a, h are the functions in (j3.27p - (j3.28p . Thus, we need to check that 

(i) IP{ ■ , 00)(x) = -b(xo,xi)Pix, 00), 

(ii) /P( • , 01)(x) = -a{xo,xi)P{x, 01), 

(iii) /P( • , 10)(x) = b{xo,xi)P{x, 00), 

(iv) IP{ ■ , ll)(x) = a{xo,xi)P{x, 01). 



(3.64) 
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Since J2yeSi ' i y) = -^1 = 0) it suffices to check only three of these equations, say (i), (ii), 
and (iv). We observe from (fOOjl - ffOT]) that /P(-,01) = -/P(-,ll). In view of this, it 
suffices to check only (i) and (ii). By I^TH}-^?^, and ([CTD - ([H:BU|) . we need 



to check that 




1 * ^ 








i * 




e 


") 


* * * ) 







/ 



(3.65) 



and 



-(1-0' -2(1-0 \ /* 1-C 2(1-0 \ /O l-i 1 



/ \** */\000 

' (3.66) 

where • denotes the componentwise product of functions and the symbol * indicates that the 
value of the functions a and h in these points is irrelevant. We see by inspection that (I3.65P 
and ()3.66p are satisfied. I 



4 Survival 

4.1 Survival bounds 

Until further notice, we continue to study the contact process on the hierarchical group Otv 
with N = 2 and its finite analogues defined in Section 13. li Our proof of Theorem [1] (b) is 
based on the following basic estimate. 

Proposition 8 (Survival bound for finite systems) Let 5 > and let {ak)k>i be non- 
negative constants. Let X^") be the ((5, qi, . . . , Un)- contact process started in Xq^^ = Sq. Then 

n-l 

p5o[j^W ^ 0] > ( J|(l - ^(A;))) e-'^('^)* (t > 0), (4.1) 

fc=0 

where 5(0) := 6, 0^(0) := Ofc (A: > 1), and we define inductively, for n > 0, 

6(n + l):=2^(n)6{n), 

(4.2) 

ak{n + 1) := ^Ok+iin) (/c > 1), 
where ^{n) := f{ai{n)/5{n)) with f as in i3.1\) . 

Proof By Lemma O one has < 0") < k all n > 0. For < ^ < i and A; > 1, let 
Pk^^ denote the probability kernel from Sk to Sk^i defined in (I331)-(l331). Let be the 

((5, ai, . . . , a„)-contact process. Applying Proposition [2] inductively, we can couple X^"-^ to 
processes 

such that take values in Sn-m, one has xj""""^ = X^'''''"'\ xf"-""^ > xj"""") 

for ah t > 0, 

p^^(n-™-l) ^ y I ^in-^) ^ ^ p„_^^(„)(x, y), (4.3) 
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and the process X^""™) is a ((5(m), ai(m), . . . , an(m))-contact process. A little thinking 
convinces us that this coupling can be done in a Markovian way, i.e., in such a way that 
(^Xin-m-i) ^ j^{n-m-i)^ is conditionally independent of 

given By this Markovian property and the definition of Pk^^{x,y), if we 

start in the initial state Xq""^ = 60, then 

m—l 

Fixf^^ =6o]=ll{l-ak)), (4.4) 

k=0 

and Xq" ™^ = with the remaining probability. Since P[X^^" = | xj."" ™^ = 0] = 1 for 
each m, we have 

P[xf ^ ^ 0] > P[xf"""^ / 0] (0 < m < n). (4.5) 

In particular, since is a Markov process in = {0, 1} that jumps from 1 to with rate 
S{n), we observe that 

n-1 

FIX^"-^ ^0]> P[xf ^ / 0] = e-'^^^^Pixi"^ = 5o] = e--^^")* " '^(^)) - (^-6) 

A;=0 

which proves (|4.1|) . I 
As an immediate corollary to Proposition [8l we obtain: 

Proposition 9 (Survival bound for infinite systems) Let 5 > and let {ak)k>i be 
nonnegative constants satisfying Y2^=i'^k < 00. Let {£,{k))k>o be defined as in Proposition\^ 
Let X he the contact process on with infection rates as in \l.b]) and recovery rate 5. Then 
the process started in Xq = 5q satisfies 

00 

P-'^iXtT^O Vt>0] > (4-7) 

A;=0 

Proof It is easy to see that the process X^") in Proposition [8] and X can be coupled such that 
X^""^ < Xt for all i > 0. Therefore ()4.7p follows from (14. ip . provided we show that 6{k) 
as A: ^ 00. In fact, it suffices to prove this under the assumption that nfclo(-'^ ~ H^)) > 0' 
otherwise (14. 7p is trivial. Indeed, nfclo(-'- ~ ?(^)) > implies that ^{k) — > as ^ 00, which 
by the fact that 

n-1 

5{n) = 5\{{2m) (4.8) 

A:=0 

implies that 5{k) — > as A: ^ 00. I 
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4.2 The critical recovery rate 

In view of Proposition [Dj we wish to find sufficient conditions for nfclo(-'- ~ ^(^)) > 0- The 
next lemma casts the inductive formula (14.21) in a more tractable form. 



Lemma 10 (Inductive formula) Let 6{n),ak{n), and ,^(n) be defined as in Propositionl^ 
and assume that the constants (afc)fc>i CLfc positive. Set e{k) := 5{k)/ai{k) (k > 0). Then 
^{k) = f{l/e{k)) and 

e{k + l) = ^^g{e{k)) {k>0), (4.9) 

where 

g{e):=4efil/e) (e > 0), (4.10) 

and f is the function defined in 

Proof It is clear from (j4.2p that ^(n) = f{l/e(n)) and 

afc(n) = 2-"afc+„ (A: > 1, n > 0). (4.11) 
Using (14. 2p once more, it follows that 



2f(n)5(n) ai(n) 6(n) , „, , . 2 "'Un+i . / n , , ^^ 

e n + 1 = = ^^4/ ai n n = ^4e n / 1/e n . 4.12 

^a2{n) a2{n)ai{n) 2 "a„+2 



The next lemma collects some elementary facts about the function g from Lemma [TOl 

Lemma 11 (The function g) The function g defined in |^.iQ[ ) is increasing on (0, oo) and 
satisfies 

g{e) = Se^ + 0{e^) as e ^ 0. (4.13) 

Proof This follows from the fact that, by Lemma El the function e ^ /(1/e) is increasing 
and satisfies 

f{l/e) =2e + 0{e^) as e ^ 0. (4.14) 



The next proposition answers the question when the infinite product in (14. 7p is positive for 5 
small enough. 

Proposition 12 (Nontrivial survival bound) Let {ak)k>o be nonnegative constants. For 
given 6 > 0, set Il{S) := nfe^o(-'^ ~ ^i^))) where the {(,{k))k>o are defined as in Proposition\^ 
Then 11(5) is nonincreasing in 6. Moreover, Il{S) > for 5 suffiently small if and only if 

oo 

2~^log(afc) > — oo for some m>0. (4-15) 

k=m 

Proof We start by showing that 11(5) is nonincreasing in S. By continuity, it suffices to prove 
this under the additional assumption that the a^s are all positive. In this case, we observe 
from Lemma [10] and the monotonicity of g that the e(A;)'s are nondecreasing in 5. Since 
(,{k) = f{l/e{k)) and / is decreasing, it follows that the ^(fc)'s are nondecreasing in 5, hence 
n(J) is nonincreasing in 6. 
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We next show that n((5) > for 5 > suffiently small if and only if (j4.15p holds. If = 
for some A; > 1, then ^{k — 1) = /(O) = ^, hence if infinitely many of the a^'s are zero then 
n((5) = for all S > 0, while (|4.15|) is obviously violated. If finitely many of the a^s are 
zero, then we may start our inductive formulas after the first m iterations, where we observe 
that 5{m) can be made arbitrarily small by choosing 5 small enough. Thus, without loss of 
generality, we may assume that the a^'s are all positive, and under this assumption we need 
to show that Il{S) > for 6 suffiently small if and only if 

oo 

J^2-'=log(afc) > -oo. (4.16) 

k=0 

It is well-known that nfclo(l " ^(^)) > if and only if ^l^o ^(^) < Using (fil^ and the 
fact that ^{k) = f{l/e{k)), it is easy to see that this is equivalent to X]fe^o^(^) ^ 
Now assume that (I4.16P holds, and, in view of ()4.13p . define {£{k))k>o by 

e(0) := — and e(k + 1) := 9^(e(k)f (k > 0). (4.17) 



ai 



Then 



_5_ 

S \2 



\ ' a4 ^ as ' ^ (TO 

More generally, it is not hard to see that 



m 

^l2) = 9f|(9e^)2(A)4^ (4.18) 



_5_\8 __ 



e{n) = (4.19) 

«n+illfc=i(afc)^ 

By Lemma [T3l below, we can choose S sufficiently small such that J2^=o^(''^) < (|4.13p . 
there exists a c > such that g{£) < 9e'^ for all e < c. By making 6 smaller if necessary, we 
can arrange that e{n) < c for all n, hence X^^o^('^) — Yl'i^=o^i^) ^ 

On the other hand, assume that Yl'^=o ^i^) < °° some S > while (|4.16p does not 
hold. Define {£{k))k>o as in (j4.17p but with the factor 9 replaced by 7. By ()4.13p . there exists 
a c > such that 76^ < g{e) for all e < c. Making 6 smaller if necessary, we can arrange that 
e(n) < c for all n, hence e(n) > e(n) for all n. Since e{n) ^ oo by Lemma [T3l below, this leads 
to a contradiction. I 



Lemma 13 (Summability) For r] > 0, set 



Fr,{n) 



"n+inLl(«fe)^ 



(n > 0). 



(4.20) 



// i4-l(^ holds, then Yl'^=o ^riin) < oo for r] sufficiently small. On the other hand, if (J^J^ 
does not hold, then lim„,^oo Pr^i^) = oo for all rj > 0. 
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Proof We start by observing that 
3m s.t. \/n > m : Frj{n) < 1 

n 

4^ 3m s.t. Vn > m : (^2" log(r/) - log(Q„+i) - ^ 2""'' log(Qfc)) < 

k=l 

3m s.t. \/n>m: I log(7?) — 2^" log(a„+i 



fc=i 



4^ log(??) -^2 Mog(afc) < 0, 



fc=i 



which is satisfied for r] sufficiently small if (j4.16p holds. In this case, we may choose ?] > 
such that K := sup„>o F^{n) < oo and observe that for any rj' < rj 



OO / 



Y^F,{n)<KY,{^) <oo. (4.22) 

n=0 n=0 ' 

On the other hand, if (|4.16p does not hold, then a calculation as in (j4.2ip shows that for all 

r] > there exists an m such that for all n > m one has Fq{n) > 1, and therefore, for any 
< ?7 < r/', 

liminf K/(n) > liminf (-) FJn) = oo. (4.23) 

n—i-oo n— >oo Tj 



4.3 Comparison argument 

Proof of Theorem [1] (b) For N = 2, Theorem [T] (b) follows from Propositions [9] and [T2j To 
generalize this to arbitrary > 2, we will use a comparison argument. 

Let N > 2 and let X be a contact process on with infection rates as in (jl.6p satisfying 

oo 

{N'y''log{ak) > -oo for some fco > 1 (4.24) 

k=ko 

where N' = N in case is a power of two and 1 < A^' < otherwise. For notational 
convenience, we set '■= ol^N^^ {k > 1), i.e., we let (7fc)fc>i denote the constants such that 
the infection rates of X are given by (compare (jl.Gp l 

= ihj G ^N, i + j)- (4.25) 

Then, by I^M), 

oo oo oo 

^ (NT' log(7fc) = (^T" log(«fc) - E (^T'k \og{N) > -oo. (4.26) 

k=ko k=ko k=ko 

We claim that we can choose n, m > 1 such that 

(^/^m < 2" < N"". (4.27) 
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If A^' = is a power of two, then this is obviously satisfied with m = 1 for some n > 1. 
Otherwise, we observe that (j4.27p is equivalent to 

log(iVO < n ^ log(iV) ^ 
log 2 m log 2 

which is satisfied for some m,n > 1 since N' < N. 

By an obvious monotone coupling, we can estimate X from below by a contact process X' 
on J7jv with infection rates of the form a'{i,j) = 7|j_j|) where 

7i = 72 = ■ ■ ■ = 7m = min{7i, • • • , 7m}, 

7m+l = 7m+2 = • • • = 72m = min{7m+i, . . . , 72m}, (4.29) 

etcetera. 

Next, we consider a contact process X" on 0,2 with infection rates of the form a"{i,j) = 7|i_j|, 
where 

// // // / 

7l = 72 = • ■ ■ = 7n = 7m, 

7n+l = 7n+2 = • • • = 72n = 72m, (4-30) 

etcetera. 

We claim that if X" survives for a certain value of the recovery rate, then so does X' . To 
see this, note that we can in a natural way identify X' with a contact process on Otv™ with 
infection rates a{i,j) = 7^|^_j|. Likewise, we can in a natural way identify X" with a contact 
process on ^2". Since 2" < A^*", we may regard r22" as a subset of Qf^fm. Therefore, by 
surpressing infections that go outside we may estimate X' from below by X" . 
For / > 0, choose ii € {Im + 1, . . . , Im + m} such that 

7i, = min{7,m+i, . . . , 7/m+m}- (4.31) 

Then, for / > and r = 1, . . . , n, one has j'/n+r ~ 7«ii hence by (14.271) and (j4.26p . 

OO OO Th OO Tl 

5] 2-'^|log(7Dl = 5]5;2-('"+'^)|log(7i;+J| = J^5^2-('"+'-)|log(70l 

k=nlo+l l=la r"=l i=Zo ^'=1 

OO OO OO 

< 5^2-'"|log(70l < Y.^N'y'"\\og{y,,)\ < (AT')'" j;(iV')"Nlog(70l (4.32) 

OO 

i=m,lo+l 

for some Iq > 0. If we write the infection rates of X" in the form a"{i,j) = aj^_^|A^^I*^'' l, then 
by and the calculation in (lOSD one has Efelfeo 

{a'D > —OO for some ko > 0, 

hence applying what we have already proved for N = 2 we conclude that X" has a positive 
critical recovery rate and the same must be true for X. I 
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A Coordinate reduction 

In this appendix we prove that Lemmas [6] and [7] imply formulas (j3.36p (i) and (ii), respectively. 
The main problem is to invent good notation. Recall that Sn = {0, 1}^". For any x £ Sn and 
A C n"", we let 

:= (x(i))ieA (A.l) 

denote the restriction of x to A. If A, A' are disjoint sets, x € {0, 1}^ and x' € {0, l}'^ , then 
we define xkx' e {0, l}^uA' 

/ / if i G A, , , 

(x dzx := < .„ . (A. 2) 

^ ^ [ x'{i) if z G A' ^ ^ 



For each i £ Q"- ^, we define Bi C by (recall ([27 

B{{) := Bi{i) = {i' oi:i' en^}. (A.3) 
Let R, R' be as in Lemma [H Then we can write 

Rf{x)= Rix,x')fix'), 



x'eSi 

R'f{y)= R'{y,y')f{y'), 



(A.4) 



where R{x,x') and R'{y,y') are the matrices of R and R', respectively. We observe that 

Rif{x)= Y Ri^\B(i),z)f{x\^-a\B(i)^z), 
.e{o,i}sw 

2e{o,i}^^> 

where we identify {0, ^ {0, 1}^' = Si and {0, l}^^ ^ {0, 1}^" = Sq. Moreover, 



Pf{x)= Y ( n p(.^^Mj)))f{y)- (A.6) 

Using the identification x\B(i) — Xi, we calculate 

RiPf{x)= Y R{xuz)Pf{x\n^\B(i)^z) 
ze{o,i}^(') 

= Y Rixi,z) Y ( n pi(^y'-\Biz)^z)j^yU)))f{y) 



Y Yl R{xi^z)piz^y{i)){ n pixj^yU)))f{y) 

^ p{xi,z)R'{z,y{i))[ n 

Y X] n pixj,yij)))fiy\n"-^\{i}^z) 

2/G{0,l}"""^ jrGO"-i 2e{0,l}{*} 

(A.7) 
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Here we have used Lemma [6] in the fourth equahty. In the fifth equahty, we have reordered 
our sums by relabelhng y{i) and z. 

The formal proof of formula (13.360 (ii) is similar, but even more cumbersome. Letting / 
and I'^ be as in Lemma [71 we can write, in matrix notation, 

If{x)= ^ I{x;x')f{x') = ^ J2 Hxo,xi;z,z')f{z,z'), 
y'eSi 2G{0,i}2'e{0,i} 



Then 



Iijf{x)= J2 E Hx\B{i),x\B(jy,Z,z')f{x\^n\(^B{i)UBij))^Z^z'), 

Ujf{x,y)= Yl E ^'x\Bu.,x\BuXy(i)^yUy,z,z')f{y\n^-i\{ij}kzkz'), 



ze{o,i}^'> z'e{o,i}{j} 



and 



(A.8) 



(A.9) 



Pfi^)= E ( n Pixk,ym)fix,y). (A.IO) 
Using the fact that — Xi, we calculate 

hjPf{x)= ^ ^ I{Xi,Xf,Z,z')Pf{x\^r.\|^B^J^yJB{j))^zkz') 

2e{o,i}s(i) z'&{o,i}BU) 
= E E I{xi,Xj;z,z') 

E ( n Pii^\^''\{B(t)vjB{j))^zkz')k,y{k))^f{y) 

= E E E Hxi,Xj;z,z')p{zi,y{i))p{zj,y{j)) 

s/e{o,i}"""^ 2G{o,i}s(0 z'g{o,i}-^(^') 

ken'^-^\{i,j} 

= E EE Pix^^z)p{xj,z')^^^.^^{z,z'■,y{^),y{j)) 

s/G{o,i}"""^ ^e{o,i}{'} ^'e{o,i}{j} (A. 11) 

n Pi^k,y{k)))f{y) 

k£n"-^\{i,j} 

= E EE pi^i^yii))pi^j^yij))^'x,,x,iyii),yU);z,z') 

•( n p(^fc>y(^)))/(yb"-i\{i,i} &^&^') 
= E ( n Pi^k,y{k))) 

■ E E lLxjiyi^)^yU);z,z')fiy\nn-i\{ij}kzSzz') 
2G{o,i}{»> 2'e{o,i}0} 
= Plijfix). 

Here we have used Lemma [7] in the fourth equality, and in the fifth equality, we have reordered 
our sums by relabelling y{i),y{j), z, and z'. 



28 



Acknowledgements 

We thank the referees for a careful reading of the manuscript, and Balint Toth, Deepak Dhar 
and Roman Kotecky for pointing out the references |Dys69| and |BM87| . 

References 

[BM87] P.M. Bleher and P. Major. Critical phenomena and universal exponents in sta- 

tistical physics. On Dyson's hierarchical model. Ann. Prohah. 15(2), 431-477, 
1987. 

[BEI92] D. Brydges, S.N. Evans and J.Z. Imbrie. Self-avoiding walk on a hierarchical 

lattice in four dimensions. Ann. Prohah. 20(1), 82-124, 1992. 

[DawOO] D.A. Dawson. Stochastic models of evolving information systems, (Ottawa, 

Canada, 1998), 1-14, CMS Conf. Proc. 26, AMS, 2000. 

[DG93] D.A. Dawson and A. Greven. Hierarchical models of interacting diffusions: mul- 

tiple time scale phenomena, phase transition and pattern of cluster-formation. 
Prohah. Theory Relat. Fields 96, 435-473, 1993. 

[DG06] D.A. Dawson and L.G. Gorostiza. Percolation in a hierarchical random graph. 

ArXiv:math/0607131v2. 

[DK96] P. Donnelly and T.G. Kurtz. A countable representation of the Fleming- Viot 

measure- valued diffusion. Ann. Prohah. 24(2), 698-742, 1996. 

[DK99] P. Donnelly and T.G. Kurtz. Genealogical processes for Fleming- Viot models 

with selection and recombination. Ann. Appl. Prohah. 9(4), 1091-1148, 1999. 

[Dur88] R. Durrett. Lecture notes on particle systems and percolation. Wadsworth & 

Brooks/Cole, Pacific Grove, 1988. 

[Dys69] F.J. Dyson. Existence of a phase transition in a one-dimensional Ising ferro- 

magnet. Commun. Math. Phys. 12, 91-107, 1969. 

[FS04] K. Fleischmann and J.M. Swart. Trimmed trees and embedded particle systems. 

Ann. Prohah. 32(3A), 2179-2221, 2004. 

[HHWOl] T. Kara, T. Hattori and H. Watanabe. Triviality of hierarchical Ising model in 
four dimensions. Commun. Math. Phys. 220(1), 13-40, 2001. 

[HL78] R.A. Holley and T.M. Liggett. The survival of the contact process. Ann. 

Probab. 6, 198-206, 1978. 

[Kur98] T.G. Kurtz. Martingale problems for conditional distributions of Markov pro- 

cesses. Electronic J. Probab. 3, Paper no. 9, 1-29, 1998. 

[Lig85] T.M. Liggett. Interacting Particle Systems. Springer- Verlag, New York, 1985. 

[Lig92] T.M. Liggett. The survival of one-dimensional contact processes in random 

environments. Ann. Prohah. 20, 696-723, 1992. 



29 



[Lig95] T.M. Liggett. Improved upper bounds for the contact process critical value. 

Ann. Probab. 23, 697-723, 1995. 

[Lig99] T.M. Liggett. Stochastic Interacting Systems: Contact, Voter and Exclusion 

Process. Springer- Verlag, Berlin, 1999. 

[RP81] L.C.G. Rogers and J.W. Pitman. Markov functions. Ann. Probab. 9(4), 573-582, 

1981. 

[SF83] S. Sawyer and J. Felsenstein. Isolation by distance in a hierarchically clustered 

population. J. Appl. Probab. 20 1-10, 1983. 

[Swa07] J.M. Swart. Extinction versus unbounded growth. Habilitation Thesis of the 

University Erlangen-Niirnberg, 2007. ArXiv:math/0702095vl. 

[Swa08] J.M. Swart. The contact process seen from a typical infected site. Pub- 

lished online in J. Theor. Probab. (2008) DOI 10.1007/sl0959-008-0184-4. 
ArXiv:math.PR/0507578v5. 



30 



